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Abstract 

There is an ever-increasing interest in investigating dynamics in time- 
varying graphs (TVGs). Nevertheless, so far, the notion of centrality in 
TVG scenarios usually refers to metrics that assess the relative importance 
of nodes along the temporal evolution of the dynamic complex network. 
For some TVG scenarios, however, more important than identifying the 
central nodes under a given node centrality definition is identifying the key 
time instants for taking certain actions. In this paper, we thus introduce 
and investigate the notion of time centrality in TVGs. Analogously to 
node centrality, time centrality evaluates the relative importance of time 
instants in dynamic complex networks. In this context, we present two 
time centrality metrics related to diffusion processes. We evaluate the two 
defined metrics using both a real-world dataset representing an in-person 
contact dynamic network and a synthetically generated randomized TVG. 
We validate the concept of time centrality showing that diffusion starting 
at the best classified time instants (i.e. the most central ones), according 
to our metrics, can perform a faster and more efficient diffusion process. 
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1 Introduction 


We witness an ever-increasing interest in investigating the dynamics of complex 
networks (i.e., changes in nodes or edges over time) representing them as time- 
varying graphs (TVGs) 11 • In this context, a number of recent efforts 


investigate new centrality definitions to capture the relative node importance in 
TVGs PPPinpp8po| . For some TVG scenarios, however, more important 
than identifying the central nodes under a given definition, is identifying the 
key time instants for taking certain actions. Deciding when, and not only where 
from, to start a diffusion process can be of upmost importance for a more effec¬ 
tive outcome. For example, Spasojevic et al. 15 has recently shown evidence 


that recommending best times for a user to post on social networks improves 
the probability of audience responses. 

This paper introduces and investigates the notion of time centrality in dy¬ 
namic complex networks. Analogously to node centrality, time centrality as¬ 
sesses the relative importance of a given time instant within a TVG, which can 
be defined in different ways. In this context, we present two metrics, namely 
cover time and time-constrained coverage, to assess time centrality in TVGs 
from different perspectives. In doing so, we also show the relevance and relative 
importance (i.e. centrality) that some time instants may reach in comparison 
with other (less important) time instants. Both metrics regard different per¬ 
spectives related to diffusion processes in complex dynamic networks. 

We evaluate both time centrality metrics we define relying on a real-world 
dataset of an in-person contact network in a hospital environment as well as on 
a synthetically generated randomized TVG (details in Section]^. Our results 
show that start spreading information at the top ranked (i.e., the most central) 
time instants, according to our metrics, can make diffusion processes in TVGs 
faster and more efficient. For example, compared with a random choice, a 
diffusion process based on time centrality metrics in the adopted real-world 
dataset can reach 10% of the studied TVG about 3 times faster in the median 
case. Likewise, a diffusion process starting at the most central time instants can 
cover more nodes in the analyzed TVGs with a limited number of allowed time 
steps for spreading when compared with random chosen time instants. 

The paper is organized as follows. Section discusses related work on tem¬ 
poral centrality in complex networks. In Section we present the model we 
use to represent TVGs. Section introduces the notion of time centrality and 
formalize two time centrality metrics related to diffusion processes. In Sectionj^ 
we discuss results on analyzing time centrality in a real-world dataset represent¬ 
ing a in-person contact network and a randomized TVG. Finally, we summarize 
results and present the outlook of future work based on time centrality in Sec¬ 
tion |6l 
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2 Related Work 


A large nuniber of works on dynamic complex networks concerns both data 
structures to model TVGs and different definitions of node centrality in that 
context. The first works that attempted to represent TVGs were influenced 
by traditional graph modeling. For instance, a number of works model TVGs 
through a sequence of static subgraphs that represents the network dynamics 
in a discrete way as time passes 


17-19 


Recent works tend to represent TVGs in different manners. For example, 
Casteigts et al. use a presence function defined over continuous time intervals. 


The model proposed by Kostakos 11 is based on the idea that a class of links 


represent instantaneous interations between distinct nodes while another class 
represents a waiting state of a given node. Kim and Anderson suggest the 
use of edges connecting nodes at different time instants (temporal and mixed 
dynamic edges). For any of these TVG models, recent proposals investigate 
adaptations of the node centrality notion from traditional graphs to the time- 
varying context, leading to different notions of temporal node centrality targeted 
at particular applications [^[^ [I2|[I4{[T^[I^[^ . 

As far as we know, this paper is the first to propose the notion of time 
centrality in dynamic complex networks. In contrast to previous work, this 
paper intends to assess the most important time instants, thus exploring the 
notion of time centrality. 


3 Modeling Time-Varying Graphs 


We model TVGs as a particular case of a MultiAspect Graph (MAG) 
in which the vertices and time instants are the key features (i.e., aspects) to 
be represented by the model. A MAG is a structure capable of representing 
multilayer and time-varying networks while also having the property of being 
isomorphic to a directed graph. The MAG structural form resembles the mul¬ 
tilayer structure recently presented by 10 , since in both cases the proposed 


structure has a construction similar to an even uniform hypergraph associated 
with an adjacency concept similar to the one of simple directed graphs. 

Formally, a MAG can be defined as an object H = {A,E), where A is a 
set of edges and A is a finite list of sets, each of which is called an aspect. In 
our case, for modeling a TVG, we have two aspects, namely vertices and time 
instants, i.e. |A| = 2. For the sake of simplicity, this 2-aspect MAG can be 
regarded as representing a TVG as an object H = (V, E, T), where V is the set 
of nodes, T is the set of time instants, and ECVxTxVxTis the set of 
edges. As a matter of notation, we denote V{H) as the set of all nodes in H, 
E{H) the set of all edges in i7, and T{H) the set of all time instants in H. 

An edge e € E{E{) is defined as an ordered quadruple e = (u, ta, v, h), where 
u,v G V{E[) are the origin and destination nodes, respectively, while ta,tb & 
T{H) are the origin and destination time instants, respectively. Therefore, 
e = {u,ta,v,tb) should be understood as a directed edge from node u at time ta 
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to node v at time tb- If one needs to represent an undirected edge in the TVG, 
both {u,ta,v,tb) and (y,tb,u,ta) should be in E[H). 

We also define a temporal node as an ordered pair {u,ta), where u G V{H) 
and ta G T{H). The set VT{H) of all temporal nodes in a TVG H is given 
by the cartesian product of the set of nodes and the set of time instants, i.e. 
VT{H) = V {H)xT{H). As a matter of notation, a temporal node is represented 
by the ordered pair that dehnes it, e.g. {u,ta)- 

The usage of the object H = {V,E,T) to represent a TVG is formally 
introduced in 23 . Therein, the representation of the TVG based on temporal 


nodes is proven to be isomorphic to a directed static graph. This is an important 
theoretical result since this allows the use of the isomorphic directed graph as 
a tool to analyze both the properties of a TVG and the behavior of dynamic 
processes over a TVG, as done is this work. This model is also shown to unify 
the representation of several previous (classes of) models for TVGs of the recent 


literature, which in general are unable to represent each other 23 


Importantly, the adopted TVG model is a particular case of a MAG 
Therefore, adding new aspects, such as different network layers, is straightfor¬ 
ward. This opens perspectives of extending the analysis of the time centrality 
concept we introduce in Section]^ to more complex dynamic networked systems, 
such as time-varying multilayer networks T pOlpT . We intend to further explore 
such perspectives in our future work. 

In the context of this paper, a time instant is equivalent to a snapshot 
in the TVG. Of course, defining the snapshot faces similar problems as any 
discretization process. For real datasets, snapshots covering large time intervals 
may not adequately capture the topology dynamics, whereas snapshots covering 
too short time intervals may generate excessive representation data without 
relevant information. To avoid these issues, in this paper, we chose to have 
snapshots corresponding to the event granularity in the case of the realworld 
dataset, i.e. 30s. Similarly, in the case of the randomized TVG, the granularity 
of the snapshots is equivalent to the changes in the topology in an atemporal 
way since we consider a new random graph in each snapshot. These TVGs used 
for the evaluation are explained in further detail in Section 


4 Time Centrality 

We here introduce the notion of time centrality. As with node centrality, dif¬ 
ferent definitions of time centrality can reflect different notions of importance 
for time instants targeting distinct applications. In other words, as node cen¬ 
trality creates a node ranking adopting an arbitrary importance criterion, time 
centrality creates a ranking of time instants using some importance criterion to 
measure the relative importance a time instant may have when compared with 
others. 

We present two time centrality metrics that can assess diffusion in TVGs 
from different perspectives. To define such metrics, we consider a generic TVG 
H = {V,E,T) (see Section]^. Let N = \T{H)\ be the size of set T{H). Each 
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time instant ti G T{H), where 0 < * < iV, composes the set of sequential time 
instants described by T[H) = For the sake of clarity, the 

illustrative figures in this section show sequences of time instants as disconnected 
snapshots, although in the adopted model temporal edges connect nodes in 
different time instants, easing the temporal analysis. 


4.1 Cover Time (CT) 

The Cover Time (CT) metric from a time centrality perspective assesses the 
minimum number of time instants (steps) a diffusion takes to reach a given frac¬ 
tion of the nodes. More formally, the cover time CT{ti,T) returns the average 
amount of time instants (steps) among all nodes for a diffusion starting at time 
ti to reach a given fraction r of the complex dynamic network, i.e. 


CT{U,t) 


{ vC7T^\ X! dt{ti,u,T), if T reached, 
' ^ '' ueviH) 

00 , otherwise. 


( 1 ) 


where dt{ti,u,T) is the number of time instants (steps) for a diffusion starting 
on node u at time instant C to reach a fraction r of nodes. If the fraction r 
of nodes is not reached for a diffusion process starting on any node at ti, then 
CT{ti,T) is considered oo. The results presented in Sectionconsider only the 
finite results of cover time. 

We consider the diffusion process to follow a Breadth-First Search (BFS), i.e. 
a generic node (u, U) starts diffusion by sending information to all of its adjacent 
nodes. Then, the adjacent temporal nodes relay information for their own adja¬ 
cent nodes in the next time instant, and so on. Information is distributed until 
the fraction t of nodes is reached and the time needed for diffusion is stored. 
This evaluation is repeated for every node {u,ti) G VT{H). The average diffu¬ 
sion time is the cover time centrality of the time instant ti. 

As an example, the top part of Figure[2shows an illustrative diffusion process 
from a single node starting at ta and finishing at ta+io, whereas the bottom part 
of Figure[2presents another diffusion process starting at ti, and finishing at tb+g. 
Intuitively, the time instant tb is more important (central) in regard of its cover 
time, since the diffusion starting at it is faster. Recall that the CT metric 
actually considers an average of the diffusion processes starting at each node. 


4.2 Time-Constrained Coverage (TCC) 

The Time-Constrained Coverage (TCC) metric from a time centrality perspec¬ 
tive evaluates the TVG coverage achieved by a diffusion process after a limited 
number of time steps. Here we assume the same diffusion process defined for the 
cover time (Section |4.1[ ). More precisely, for a given time instant U, TCC{ti, $) 
returns the average fraction of the TVG nodes reached by the diffusion process 
in $ of time steps. More formally, 
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Figure 1: Illustrative cover time on TVGs. 


TCC{U,^) = —^ ^ ( 2 ) 

where dc(ti, u, <&) is the number of nodes reached from node u after $ time steps 
when the diffusion starts at time ti. For example, Figure [^presents two distinct 
illustrative diffusion processes occurring in $ time steps from a single node. 
At the left part, the process starting at ta reaches less nodes than the process 
starting at at the right part after $ time steps. The latter diffusion process 
actually reaches all nodes in the TVG at the time instant. Therefore, the 
time instant tb is more important (central) in regard of its time-constrained 
coverage, since the diffusion starting at it is more effective in reaching a larger 
number of nodes after $ time steps. Again, recall that the TGG metric considers 
an average of the diffusion processes starting at each node. 


5 Time Centrality Results 

In this section, we first present the TVGs we use to evaluate the two time 
centrality metrics. We next analyze the results for the cover time metric. Then 
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Figure 2: Illustrative time-constrained coverage on TVGs. 


we analyze the results for the time-constrained coverage metric. Finally, we 
further discuss on the achieved outcome of the time centrality metrics as well 
as on some limitations of the study. 

5.1 Real-world dataset and randomized TVG 

We evaluate the time centrality metrics we propose using a TVG dataset col¬ 
lected in the context of the MOSAR project [^. MOSAR (Mastering hOSpital 
Antimicrobial Resistance and its spread) is a scientific collaboration project 
that comprises several medical, biochemistry, and computing research institu¬ 
tions. The MOSAR project focuses on antimicrobial-resistant bacteria (AMRB) 
transmission dynamics in high-risk environments, such as intensive care units 
and surgical centres. The adopted dataset consists of the records of in-person 
contacts (from physicians, nurses, staff members, and patients) in a certain med¬ 
ical ward for a period of two weeks (between 12 am of July 25, 2009 and 12am of 
August 08, 2009). Each one of the 160 volunteers who participated in the study 
was equipped with a RFID device that detected the presence of another devices 
within a small distance (about one meter). Device identification was unique 
and was always associated with the same person. Every 30 seconds during the 
two-weeks period, each device registered the list of all devices (nodes) that were 
within its coverage area in order to establish the arrangement of the contacts 
among them (edges). 

We have modeled the MOSAR dataset as a TVG M = (V, E, T), as described 
in Section]^ where V{M) is the set of all participants, E{M) is the set of all 
in-person contacts, and T{M) is the set of sequential time instants at a 30 
seconds granularity. Thus, |V(M)| = 160 and \T{M)\ = 40320. Gontacts are 
represented by non-oriented edges, i.e., the edge between any two persons is 
represented by a directed edge and its reciprocal. 

We also evaluate the time centrality metrics using a randomized TVG R = 
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{V,E,T), where V{R) is the set of nodes, E{R) is the set of edges, and T{R) 
is the set of time instants, represented by a sequence of independent snapshots. 
Each snapshot is synthetically generated as a random graph based on the Erdos- 
Renyi G = (n^p) model, where n is the number of nodes in the graph and 
p is the probability for edge creation between each pair of nodes. For the 
randomized TVG, we have kept the same number of nodes of MOSAR TVG, i.e., 
|E(i?)| = n = 160, and we have arbitrarily defined the number of time instants 
jT(i?)| = 800. Furthermore, we have defined the probability p = 0.01 x 
ensuring that each resulting random graph corresponding to a snapshot presents 
a similar level of sparsity as compared with the MOSAR TVG (i.e. the value 
0.01 ensures this sparsity). This value of probability p also ensures that the 
generated random graphs are disconnected since it is known that p > is a 
sharp threshold for the connectedness of G = (n.,p) random networks [^. 


5.2 Analysis of Cover Time 

We first present results for the cover time metric. As previously discussed, lower 
cover time values indicate more central time instants that can spread informa¬ 
tion in a shorter period of time. Figure 3(a) and |3(b)] present the cumulative 
distribution function (GDF) of the cover time for different values of the frac¬ 
tion T of nodes to be covered in the MOSAR TVG and in the randomized 
TVG, respectively. Six different scenarios in each subfigure are considered with 
T = {0.1, 0.2,..., 0.6}. In Figure |3(a)| results compare the first 33000 time 
instants, so there is room for the diffusion to spread through the 40320 time 
instants of the studied MOSAR TVG. For the reason, in Figure |3(b)[ results 
compare the first 660 time instants of the 800 time instants composing the 
randomized TVG. 

In Figure]^ each value in a curve represents the fraction of the studied time 
instants that achieved the corresponding t fraction of nodes in at most the 
corresponding number of time steps, i.e. the cover time t. In order to improve 
visualization quality, we add symbols on arbitrary positions of the curves (the 
same applies to all remaining figures). For instance, in Figure [3(a)[ to cover 
60% of the MOSAR TVG (r = 0.6), 20% of the time instants take less than 
2000 time steps. In other words, information diffusion starting at 20% of the 
time instants takes less than 16 hours and 40 minutes to cover 60% of the nodes 
in the MOSAR TVG only spreading opportunistically through the in-person 
contact dynamic network. In Figure [3(b)[ to cover the same fraction of nodes, 
20% of the time instants take less than 100 time steps. 

From Figure]^ for both cases, we clearly note that starting the diffusion at 
a few central time instants reduces significantly the time of the spreading to 
reach a certain portion of the nodes in the TVG. Therefore, we also measure 
the relative improvement of diffusion efficiency achieved by using only the most 
central time instants. In Figure]^ we compare for both cases, i.e. the MOSAR 
TVG and the randomized TVG, the performance of the top 10 most central time 
instants according to the cover time metric against 10 randomly chosen time 
instants to start spreading a diffusion. Taking the mean result as a reference, to 















Cover Time t (# of time instants) 


(a) Cover time distributions in the MOSAR TVG for varying values 
of the fraction r of nodes to be covered. According to the results, to 
cover 10% of the network (r = 0.1), the median cover time is 513 time 
steps. Similarly, to reach 60% of the network (r = 0.6), the median 
cover time is 2994 time steps. 



(b) Cover time distribution in the randomized TVG for varying values 
of the fraction r of nodes to be covered. According to the results, 
to cover 10% of network (r = 0.1), the median cover time is 50 time 
steps. Similarly, to reach 60% of the network (r = 0.6), the median 
cover time is 106 time steps. 

Figure 3: Cover time distributions in MOSAR and randomized TVGs. 
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(a) Cover time in the MOSAR TVG. To reach 10% of the network (r = 
0.1), the median cover time is 206 time steps, if the diffusion starts at 
the top 10 time instants; whereas the median cover time is 815 time 
steps, if the diffusion starts at 10 randomly chosen time instants. 



(b) Cover time in the randomized TVG. To reach 10% of the net¬ 
work (r = 0.1), the median cover time is 27 time steps, if the diffusion 
starts at the top 10 time instants; whereas the median cover time is 48 
time steps, if the diffusion starts at 10 randomly chosen time instants. 

Figure 4: Cover time in MOSAR and randomized TVGs: top 10 vs Random 
Choice. 
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cover 10% of the nodes (t = 0.1) in the MOSAR TVG in Figure [4(a)| a diffusion 
starting at the top 10 most central time instants takes 3 times less time steps 
than a diffusion starting at the random set of time instants. Similarly, for the 
randomized TVG in Figure [4(b)[ a diffusion starting at the top 10 most central 
time instants takes almost half the time steps than a diffusion starting at the 
random set of time instants to cover 10% of the nodes (r = 0.1). 

5.3 Analysis of Time-Constrained Coverage 

The time-constrained coverage is also able to select few time instants that 
present distinct diffusion performance. Figure presents the complementary 
CDF of the time-constrained coverage of the MOSAR and randomized TVGs 
for different values of $ time steps. Likewise in cover time of the MOSAR TVG, 
the first 33000 time instants are compared on each case, while the first 660 time 
instants are compared on each case of the randomized TVG. 

In Figure each value in a curve represents the fraction of the studied time 
instants that achieved at least the corresponding time-constrained coverage c of 
the TVG in $ time steps. For example, in Figure |5(a)[ with a limit of <i> = 2000 
time steps diffusion starting at 80% of the time instants reach at least 28% of 
the nodes in the MOSAR TVG, noting that the maximum reachable coverage 
is 58%. 

As expected, the larger is the maximum number $ of allowed time steps, 
the larger is the overall network coverage. Importantly, however, in all cases 
a very limited number of time instants — actually, the most central ones — 
provide superior efficiency in covering more nodes in both the MOSAR and the 
randomized TVGs with a constrained number of time steps when diffusion starts 
at these most central time instants. In Figure for both the MOSAR and the 
randomized TVGs, we compare the performance of the top 10 most central time 
instants according to the time-constrained coverage metric against 10 randomly 
chosen time instants to start spreading a diffusion. Indeed, in Figure |6(a)[ 
starting a diffusion on any of the top 10 most central time instants would reach 
at least 56% of network coverage, whereas in the remaining randomly chosen 
time instants, a diffusion process would reach no more than 33% of the MOSAR 
TVG nodes, considering a limit of <I> = 2000 time steps. 

5.4 Further Discussion 

Overall, the observed time centrality results, for all considered scenarios, do 
not follow a normal or any other narrow distribution. Indeed, this has been 
evidenced in Figures]^ andwhere we clearly note that random choices of time 
instants present results for the considered metrics that are significantly worse 
than the top 10 ranked time instants. This kind of behavior typically indicates 
heavy-tailed distributions, thus highlighting the relevance of time centrality in 
distinguishing the relative importance of different time instants for diffusion 
purposes. 
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In the experiments, both the MOSAR and the randomized TVGs have the 
same number of nodes and a similar node density. Despite that, we have ob¬ 
served that the time centrality results differ in magnitude when we compare 
both the results for the MOSAR and the randomized TVGs, although they 
present a similar general behavior. This difference in the magnitude of results is 
a consequence of the contact opportunities nodes have in the considered TVGs. 
In fact, in the randomized TVG, more than 99% of edges shift from connected 
to disconnected (and vice-versa), thus providing uniform contact opportunities 
for virtually all nodes, allowing faster as well broader diffusion. In contrast, the 
dynamic in-person contact network represented by the MOSAR TVG presents 
only 56% of edges shifting their current states, thus significantly limiting the 
contact opportunities nodes have to diffuse information. As a consequence, dif¬ 
fusion becomes much slower and constrained in the MOSAR TVG as compared 
with the randomized TVG. 

Finally, we highlight that, in this paper, the intended focus was mainly on 
pointing out the relevance of the time centrality concept. We thus expect this 
paper to motivate further research on how to detect central time instants us¬ 
ing past or present information. This opens perspectives to investigate suitable 
metrics for such an evaluation as well as detection systems of central time in¬ 
stants based on those metrics. These natural next steps in the way the paper 
points out are our intended future work. 


6 Summary and Outlook 

In this work, we introduce the notion of time centrality in dynamic complex 
networks. Time centrality assesses the relative importance of a given time in¬ 
stant within a time-varying graph (TVG). In this context, we present two time 
centrality metrics focused on diffusion processes and evaluate them using a real- 
world dataset representing a dynamic in-person contact network. Our results 
show that starting a diffusion at the most central time instants, according to 
our centrality metrics, can make diffusion processes in TVGs faster and more 
efficient. 

Considering the notion of time centrality in dynamic complex networks opens 
several perspectives for further research. As future work, we plan to further 
analyze the particular case of considering time centrality conditioned upon the 
seed node that holds the information to be diffused. Furthermore, we also intend 
to evaluate the trade-off on the diffusion performance in TVGs between waiting 
for the most central time instant to spread information or starting the diffusion 
before the best moment with better performance expectations considering the 
total time of both retaining and spreading of information. 

Moreover, in our future work we plan to develop prediction models based 
on time centrality for complex systems that can be represented by TVGs. In 
that sense, strategies for the possible early identification of central time instants 
might be based on discovering evidence in the recent past or the present moment 
of the TVG evolution indicating that a given time instant is relatively better 
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(i.e. more central) with respect to a set of others. 
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(a) Time-constrained coverage distribution in MOSAR TVG. Accord¬ 
ing to the results, for a $ = 500 time steps diffusion, the median 
time-constrained coverage is 10% of network. Similarly, for $ = 2000 
time steps, the median fraction of network reached is 41.5%. 



(b) Time-constrained coverage distribution in randomized TVG. Ac¬ 
cording to the results, for a diffusion of 4> = 25 time steps, the me¬ 
dian time-constrained coverage is 3.1% of the network. Similarly, for 
$ = 100 time steps, the median fraction of network reached is 43.7%. 

Figure 5: Time-constrained coverage distributions in MOSAR and randomized 
TVGs. 
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(a) Time-constrained coverage in the MOSAR TVG. The median time- 
constrained coverage for a <1> = 2000 time steps diffusion is 56%, if the 
diffusion starts at the top 10 time instants; whereas for 10 randomly 
chosen time instants, time-constrained coverage is 36%. 
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(b) Time-constrained coverage in the randomized TVG. The median 
time-constrained coverage for a 4> = 100 time steps diffusion is 60%, if 
the diffusion starts at the top 10 time instants; whereas for 10 randomly 
chosen time instants, time-constrained coverage is 52%. 

Figure 6: Time-constrained coverage in MOSAR and randomized TVGs: top 
10 vs Random Choice. 
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